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ABSTRACT 

The Wolf-Rayet star WR 46 is known to exhibit a very complex variability pattern 
on relatively short time scales of a few hours. Periodic but intermittent radial velocity 
shifts of optical lines as well as multiple photometric periods have been found in the past. 
Non-radial pulsations, rapid rotational modulation or the presence of a putative low- 
mass companion have been proposed to explain the short-term behaviour. In an effort 
to unveil its true nature, we observed WR 46 with FUSE (Far Ultraviolet Spectroscopic 
Explorer) over several short-term variability cycles. We found significant variations on 
a time scale of ~8 hours in the far-ultraviolet (FUV) continuum, in the blue edge of 
the absorption trough of the O VI AA1032, 1038 doublet P Cygni profile and in the 
S VI AA933, 944 P Cygni absorption profile. We complemented these observations with 
X-ray and UV light-curves and an X-ray spectrum from archival XMM-Newton (X-ray 
Multi-Mirror Mission - Newton Space Telescope) data. The X-ray and UV light-curves 
show variations on a time scale similar to the variability found in the FUV. We discuss 
our results in the context of the different scenarios suggested to explain the short-term 
variability of this object and reiterate that non-radial pulsations is the most likely to 
occur. 

Subject headings: stars: Wolf-Rayet - stars: individual (WR 46) - ultraviolet: stars - 
X-rays: stars 
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1. Introduction 

The Wolf-Rayet (WR) star WR 46 (HD 104994, DI Cru) is well-known for its remarkable short- 
term variability. All the suspected causes of this variability (a close, low-mass binary companion, 
rapid rotational modulation, non-radial pulsations) represent very rarely observed phenomena for 
Wolf-Rayet stars. The confirmation of either of these scenarios would therefore provide insight into 
our understanding of the physics and evolution of massive stars. 

Despite the fact that some scenarios have been favoured to explain the puzzling behaviour of 
this star, we still lack a satisfying explanation of the underlying mechanism. Because of the detec- 
tion of radial velocity variations, the presence of an unseen close binary companion was initially a 
popular suggestion. Although it runs into many difficulties, this possibility has never been com- 
pletely excluded up to now. It has also been suggested that mult i-frequency non-radial pulsations 



of a single WR star could cause the observed wind variability (jVeen et al.l l2002d ; lOliveira et al 



2004). While this last scenario appears to be at least part of the solution, more observations in 



selected spectral windows were needed to progress towards a clearer picture of this enigmatic star. 

It is for this reason that we obtained a time-resolved series of FUSE (Far Ultraviolet Spec- 
troscopic Explorer) observations and that we retrieved from the archive available XMM-Newton 
(X-ray Multi-Mirror Mission - Newton Space Telescope) data of WR 46. FUSE observed in the 
far-ultraviolet (FUV) spectral range, giving access to P Cygni profiles of resonant transitions from 
various species. These profiles are characteristic of dense stellar outflows and the study of their 
variability enables us to probe the evolution of the wind structure of WR 46 as a function of time. 
We were particularly interested in the high ionization potential O VI AA1032, 1038 doublet, which, 
along with the X-ray data from XMM-Newton, could provide information on energetic phenomena 
such as accretion onto a compact companion or wind-embedded shocks. 

Before we present these observations, we first highlight in $2] key aspects of the rich but complex 
observational history of WR 46 in order to place our own data set in the general context. In $3] we 
describe our observations and the data reduction procedures. In $4j we present our analysis of the 
data, and in $5] we discuss our results. Our conclusions can be found in $6l 



2. Observational history 



2.1. Physical properties 



WR 46 is classified as WN3p in the Seventh Catalogue of Galactic WR Stars (jvan der Hucht 



2001 



). It is characterized by the presence of strong N V and He II lines and the absence of hydrogen. 
The "p" stand s for peculiar and den otes the presence of unusually strong O VI AA3811, 3834 emis- 
sion lines (e.g. ICrowther et al.lll995T). who se occurence is probably due to the susceptibility of that 
doublet to low density (|Veen et al.1 l2002bT ) . Althoug h it has once been con s idered a low mass X-ra y 
binary ( Niemela et al. 19951 ) or a V Sagittae star ( Steiner Diaz 1998 ). Crowther et al. ( 1995 ). 



Marchenko et al. 
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Veen et al. 


(2002bllcl). 


Oliveira et al. 


(2004 


) and 


Gosset et al. 


(2011) 
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convincing ar guments in favour of the Wolf-Rayet nature of WR 46. As a weak-lined WN early 



(WNE) star (jCrowther et al.l Il995l ). it exhibits a low degree of atmospheric extension (i?* 



R TR=2 /^) and triangular emission-line profiles indicating very low wind densities. Weak-lined WNE 
stars like WR 46 are as hot as the strong-lined WNE stars, but with mass-loss rates about an order 
of magnitude lower. 



WR 46 is a known X-ray source whose dis covery with Einstein was reported by Polloc k! (119871 ) . 
Later, based on ROSAT PSPC observations ( Wessolowski et al. 1995 ). Crowther et al. ( 19951 ) es- 
timated the X-ray lumin osity of WR 46 to be Lx = 5±2x 10 31 erg s _1 f or a distance of 4.0 kpc 
(|Tovmassian et al.lll996l ) and noted that this is typi cal of single WN sta rs (|Poilocklll995l ). In their 
detailed analysis of XMM-Newton data of WR 46, I Gosset et al. ( 2011 ) deduced a luminosity of 



L x (0.2-10.0 keV)= 7.7 x 10 32 erg s" 1 (for d =4.0 kpc). Using our own X-ray model (see §4.4jl 
shows the previous Einstein and ROSAT X-ray measuremen t s to b e consistent with the X-ray 
luminosity from XMM data within about 10%. iGosset et al.1 (|201ll ) interpreted the X-ray spec- 
trum of WR 46 as dominated by a soft component but also detected a h arder component abo ve 



3 keV, similar to t he X-ray spectrum of other presumably single WN stars (jSkinner et al. 
Ignace et al.ll2003l ). 



2002 



In the present work, we adopt the results from a recent CMFGEN (IHillier Sz Millerlll998l ) blan- 
keted wind model of WR 46 (Crowther, private communication). The stellar and wind parameters, 
estimated from a combination of nitrogen and oxygen optical/UV line diagnostics, are presented 
in Table [TJ Note that for the wind terminal velocity, we adopt the value of 2775 km s _1 measured 
from the saturated O VI AA1032, 1038 P Cygni absorption line (jWillis et al.ll2004l ). a resonance 
doublet that in principle allows to trace material further out in the wind (at lower densities, closer 
to the actual terminal velocity) than the non-resonance optical/ UV lines used in othe r stud ies. The 
adopted mass was derived using the mass-luminosity relation of lSchaerer &: Maederl (| 19921 ). 



2.2. Variability 



Since iMonderen et al.l (|1988l ) first reported the photometric variability of WR 46, many studies 
have presented time-resolved photomet ric and/or spectroscopic follow-ups of this star. We refer 



to some of the more recent studies (e.g. IVeen et al.ll2002al : lOliveira et al 



20041 : IGosset et al.1 12011 



for a detailed review of the variability observed in WR 46. Here, we just outline the main periods 
found by various authors and the conclusions that were drawn about the origin of the variability. 
Figure Q] shows the many photometric and spectroscopic (radial-velocity) periods reported, along 
with the time scales determined in this work (see §4.31 and £ j4.5p . The spectroscopic "period" from 
this work is not a radial-velocity period, but it is instead the time scale of variations found in the 
FUV O VI AA1032, 1038 doublet P Cygni profile. The different periods and time scales observed 
suggest that the period of WR 46 is changing and/or that multiple periods are present. However, 
we stress that, ideally, the reality of each of the secondary (i.e. non-dominant) periods shown on 
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Figure Q] would need to be confirmed with a long-term and very intense data set. 



Although the first observations of short-term variability in WR 46 (binary-like double-wave 
light-curve, low-amplitude radial- velocity variations with K ^100 km s" 1 ) led to the suggestio n 



that it was a binary with a low- mass companion (e.g lvan Genderen et alJll99ll ; IVeen et al.lll995l ). 
subsequent observations reall y casted doubt on this hypothesis. The fact that the dominant period 
appears to change over time (jVeen et al.ll2002al ) or that the radial-velocity variations unexpectedly 
disappear in about one out of three consecutive nights (jMarchenko et al.ll2000l ; I Veen et al.ll20 02b) 
is difficult to reconcile with a binary scenario. While not complete l y discarding the po ssibility of 
a close binary, the most recent studies of WR 46 (|Veen et al.ll2002a ; lOliveira et al.ll2004l ) proposed 
that non-radial pulsations may well play a dominant role in the short-term variability of this star. 

In addition to its short-term cyclical variability, WR 46 is known to exhibit significant vari- 
ability on a time scale of months. A brightness increase occurred between 1989 and 1991, followed 
by a dec rease from 1 991 to 1993 (jMarchenko et alJll99a : IVeen et al.l l2002al ) . From a series of IUE 



spectra, IVeen et al.l (|2002bl ) in ferred that the 1991 brightness maximum was accompanied by an 



increase of the mass-loss rate. lOliveira et al.l (|2004l ) found evidence for a difference in the degree 
of ionization of WR 46 between June 1996 and January 2002 and suggested that this is probably 
caused by a variation of the global density of the wind. The link between the possible evolution of 
the short-term period(s) and the long-term brightness and global wind density variations however 
remains unclear. 



Observations and data reduction 
3.1. FUSE 



The FUSE scientific instrument has been described in detail by iMoos et al.l (120001 ) . The in- 
strument imaged four spectra (SiCl, SiC2, LiFl, LiF2) split in two segments (A and B), for a 
total of eight segments covering the spectral window from 905 to 1187 A with substantial overlap 
between individual segments. 

We secured FUSE observations of WR 46 under guest observer time (program E113, PI: N. 
St-Louis) in order to carry-out time-resolved FUV spectroscopy over many cycles of variability 
of the star. In total, for the three main consecutive observations of this program in March 2006 
(E1130103, E1130104, E1130105), the star was followed for more than 30 h ours spread over a 



timespan of about 50 hours. A previous observation of WR 46 was presented by I Willis et al.l ((2004) 
in their FUSE atlas of WR stars, and a few other shorter observations of this star are also found 
in the FUSE archive. We retrieved all available science observations of WR 46 for the analysis 
presented in this work. The log of the FUSE observations is presented in Table El which gives 
the FUSE data set number, the date and time of the beginning and end of the observation, the 
aperture and mode used, as well as the number of exposures and the total exposure time of each 
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observation. Note that the observations were obtained in time-tag (TTAG) mode, for which the 
arrival times of photons are recorded. 

Each individual ra w exposure was ru n through the latest version of the FUSE calibration 



pipeline, CalFUSE v3.2 (IDixon et al.1 120071 ). The pipeline first corrects the input data for "gaing- 
sag", "event bursts", and positional shifts of detected photons due to various effects. Then, a 
correction for background light is applied and a one-dimensional spectrum is extracted. Finally, 
the eight channel spectra are calibrated in wavelength and flux. In order to maximize the time 
coverage, we did not mask the data obtained during the daytime portion of the orbit as is often 
done for FUSE data. Geocoronal airglow emission lines are stronger for day-time data, but this 
did not represent a serious constraint here. The strongest airglow line, H I Lyman-/3 at 1026 A, is 
located right in the saturated part of the absorption component of the O VI AA1032, 1038 P Cygni 
profile and therefore it does not mask any precious spectral information. The other airglow lines 
seen in our data were much weaker. 

Since we wanted photometric exposures to later extract accurate light-curves, we modified 
CalFUSE output intermediate data files (IDF) to remove the time intervals where the count rate 
in the spectral extraction window drops relatively abruptly due to the drifting of the star out of 
the aperture or away from the aperture center. It is well known that the photometric accuracy 
of FUSE is influen ced by various effects that cannot be fully corrected by the CalFUSE pipeline 



(jDixon et al.ll2007l ). A target centered in an aperture of the guide channel (LiF 1 prior to July 2005, 
LiF 2 after that) can be misaligned in the apertures of the other three channels. Also, with the loss 
of reaction wheels (the first two in 2001), channel drifts can even temporarily move the target out 
of the guide-channel aperture, as happened for some exposures of WR 46. The pipeline attempts 
to flag times when the target is out of the aperture, but this time lost to pointing excursions is by 
default underestimated to avoid the rejection of good data. A careful examination and screening of 
the IDF allowed us to insure the photometric quality of the data. New spectra were then extracted 
from these IDF for each exposure separately. The individual exposure times of these spectra are 
typically on the order of a thousand seconds. We attempted to further break down all the exposures 
in small pieces of a few hundred seconds, but came to the conclusion that the marginal improvement 
in time resolution was generally not worth the associated decrease in the signal-to-noise ratio of 
spectra. 

The spectra were extracted from the same modified IDF that were used to obtain light-curves. 
Because the count rate during the masked intervals is lower than normal and because these masked 
intervals usually represent a minor fraction of the exposure time (especially for the guide channel 
apertures), the fraction of spectral photons lost due to the screening of the IDF is small. Thus, 
the signal-to-noise ratio of the extracted spectra is not significantly affected. In a few cases, a 
major fraction if not all of the exposure had to be masked and no spectrum was extracted, but 
then the original spectrum that was extracted previous to screening for pointing losses had a very 
low signal-to-noise ratio and was hardly usable anyway. Using photometric data from the modified 
IDF to extract spectra advantageously gives a more reliable absolute flux calibration, thus making 
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it possible to study long-term changes in the FUV spectrum of WR 46. 

When possible, an instrumental artifact nicknamed the "worm" was removed from the data. 
This feature is caused by electron repeller grid wires that prevent light from the target from reaching 
the detector, resulting in a wide artificial absorption trough. It was most noticeable in the LiFIB 
spectral segment of the observations we retrieved. It was also sometimes seen in the LiF2A segment. 
To remove the LiFIB (LiF2A) worm of a given exposure while keeping intact the smaller scale 
spectral features of the LiFIB (LiF2A) spectrum, we multiplied the LiFIB (LiF2A) spectrum by 
the smoothed ratio of the overlapping LiF2A (LiFIB) and LiFIB (LiF2A) spectra. This method 
was not always possible, as for some exposures a worm was apparent in both LiF2A and LiFIB. In 
these cases, the LiFIB and LiF2A spectral segments of these exposures were not considered in our 
analysis of variability. 

To first isolate spectral line profile variability from continuum variability, all spectra were 
scaled to a common continuum level. The spectra were also co-aligned to the spectrum of a refer- 
ence exposure by cross-correlating them over small spectral ranges that include narrow interstellar 
absorption lines. The alignment required for the different spectra was typically of the order of 0.04 
A or less. All the spectral segments from the reference exposure were previously co-aligned on the 
guide-channel (LiF 2) spectral segments, for which the wavelength calibration errors are minimal 
(jDixon et al.l 120071 ). For each exposure, we then performed a weighted average merging of the flux 
data in all eight channels, with the weights inversely proportional to the square of the statistical 
errors. In the end, we obtained a single 1-D, heliocentric corrected, flux-calibrated spectrum for 
each exposure (except when the exposure suffered from serious misalignment of the target). 

Light-curves were obtained to investigate the FUV continuum variability. Other light-curves 
were also extracted to look at variable spectral-line regions (see £ j4,2p integrated over broader 
wavelength intervals than the typical resolution of FUSE spectra, but with a much better time- 
resolution than it would be possible by just looking at our time-series of spectra. The routine 
ttagJightcurve (see Bernard Godard's IDF cookbooHj) was used to extract light-curves from the 
modified IDF discussed above. We estimated 1-a uncertainties on the data points of the light-curves 
by assuming Poisson statistics (a = y/N, where ./V is the number of counts). 

Specifically, three different types of light-curves were produced. The first shows the variations 
in the FUV continuum. The spectral window used for these FUV continuum light-curves includes all 
wavelength intervals that do not encompass the e mission lines and P Cy gni profiles identified in §4.1^ 



the strongest predicted airglow emission lines (see lFeldman et al.ll200ll ). the edges of t he wavelength 



range of each spectral segment, where systematic errors in the flux are significant (jSahnow et al 



2000), and narrow interstellar absorption features. The two other types of light-curves show the 



variations in the O VI AA1032, 1038 and S VI AA933, 944 doublet P Cygni absorption components, 
the only two variable spectral regions in the FUSE spectrum of WR 46 (see §4.21 ) . The O VI 



x http: //archive. stsci.edu/fuse/analysis/idfcook/ 
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light-curve was extracted using the wavelength region from 1019.5 to 1022.5 A, and the S VI light- 
curve the wavelength regions from 924.5 to 929 A and 934 to 939 A. These light -curves are slightly 
affected by continuum variability, but they are largely dominated by the line variability, which is 
much stronger than the continuum contribution in these selected spectral regions. This was easily 
checked, for example by realizing that the variability in these spectral regions was roughly the same 
for the spectra that were scaled to a common continuum level compared to the unsealed series of 
spectra. 

For each extracted light-curve, the contributions from all possible spectral segments were 
combined. For the FUV continuum light-curve, we combined the continuum light-curves of all 
eight FUSE spectral segments since continuum variations appeared to be the same across the 
entire FUSE spectral range. For the O VI AA1032, 1038 P Cygni absorption light-curve, all the 
spectral segments covering the spectral range concerned (namely LiFIA, LiF2B, SiCIA, and SiC2B) 
were used. Similarly, spectral segments SiClB and SiC2A, covering the absorption component of 
the S VI AA933, 944 P Cygni profile, were used for the third light-curve. Note that when a worm 
feature was present in a spectrum, the light-curve of the affected spectral segment and exposure 
was not used. 

To combine the light-curves of individual spectral segments we first normalized their average 
count rate to 1 for each observation and rebinned them into time intervals of 0.2 hours, propagating 
Poisson l-o" errors accordingly. An average light-curve of all the individual spectral segments 
concerned, weighted using errors, was then computed. To double-check the photometric quality of 
the data, we made sure that the light-curves from individual channels appeared correlated before 
averaging them. 



3.2. X MM- Newton 



WR 46 was observed by XMM-Newton (jJansen et al 



2001 



Lj) between 02:05:52 and 23:22:57 on 
2002 February 8 (ObsID 0397-0109110101; 2XMM catalogue source name = 2XMM J120518.7- 
620310) with the EPIC and OM instruments. As d iscussed in fTL| the XMM-Newton X-ray 
spectrum of WR 46 is presented bv lGosset et al.l (|201ll ). along with X-ray light-curves in different 
bands extracted from the XMM-Newton data. To inform our discussion of the short-term variability 
of WR 46, we also retrieved the publicly available XMM-Newton data of this star, and performed 
the additional task of extracting UV light-curves simultaneous to the X-ray light-curve. 

The EPIC o bservations presente d here were obtained with the pn detector used in the extended 
full frame mode (Striider et al.ll200ll ) and with the two MOS detectors used in the full frame mode 



( Turner et al 



2001 



). A medium filter to reject parasite optical/UV light from the target was used 
for all three detectors. The OM instrument was used in the imaging mode with the broadband 
ultraviolet filters UVM2 and UVW2 (effective wavelength 231 and 212 nm respectively). 



Custom X-ray source and background spectra were extracted from the calibrated event lists 
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with standard procedures in the XMM-Newton Science Analysis System (SAS) v9.0.0, using circular 
source and annular background selection regions. The light-curve was calculated with a PSF 
modeling procedure to estimate simultaneously source and background variations, whose results 
agreed well with the more approximate standard methods. This light-curve includes the whole 
XMM-Newton energy band. Although we decided to calculate our own spectra and light-curves, 
they proved to be consistent with those available through the XMM-Newton Science Archive via 
the 2XMM catalogue. 

The OM observations consisted of 25 exposures with the UVM2 filter and 30 exposures with the 
UVW2 in image mode with 1000 s exposure time for all of them except the last 10 exposures with 
UVW2 that had only 860 sx 5 and 800 sx 5. Aperture photometry of the target was performed using 
the SAS. The source was automatically detected by the SAS in all images and the count rate was 
measured using an aperture of radius 35 pixels (circa 17.5 arc sec). The background was measured 
in an annulus centered on the target. The measured rates were corrected from coincidence losses 
and time sensitivity degradation of the detector by the SAS. 

In order to later analyze the full OM UV light-curve, we put one filter into the system of the 
other using the procedure described below. To perform this correction, we have to assume that 
the total flux ratio in the band pass of the UVM2 and UVW2 filters is constant throughout the 
observations. This assumpt ion seems reasona ble given that the color variations in the VBLUW 



photometric monitoring of (jVeen et al.ll2002al ) are smaller than 1%, and the a mplitude of photo 



metric variability is around ±5% in both the UV (see §4.5h and the optical (e.g. lVeen et al.ll2002al ). 
We then have: 

-FUVM2 

where F refers to the flux and K is a constant. From the calibration of the OM, we know that: 



^filter — -Rfilter Cfilterj (2) 

where F is the flux, R is the count rate and C is a rate to flux conversion factor. This conversion 
factor can be obtained from the effective area and observations of standard stars, from the OM 
calibration documentation, or from the SAS Current Calibration File (CCF). 

To compute K, we used the series of observations of WR 46 in the UVW2 and UVM2 filters 
as follows: 

K _ Ryvw2 Cuvw2 

-RtJVM2 CtjVM2 

where -RtjvW2 an d -Rtjvm2 are the mean count rates in each filter. The normalization to convert 
the measured rate in the UVM2 filter to the UVW2 equivalent is then 
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-RUVW2 = K -RlJVM2 UVM2 , (4) 
<^UVW2 

Here, we have -Ruv\V2 = 152.417 counts s -1 , -RuvM2 = 292.748 counts s _1 , Cuvw2 = 5.75 x 
10~ 15 erg/cm 2 /A/count, and CuvM2 = 2.20 x 10~ 15 erg/cm 2 /A/count. The same method was 
applied to compute normalized errors. 



4. Analysis 
4.1. Time-averaged FUSE spectrum 



From a visual inspection of the time-averaged FUSE spectrum of WR 46 (see Figure [21 top 
panels) an d the synthe t ic int erstellar transmission spectra for the FUSE range computed, for ex- 
ample, by IWillis et al.l ([2004), it is clear that the spectrum of WR 46 is highly contaminated by 
ISM features. The column density of interstellar H2 is perhaps as high as 10 18 atoms/cm 2 and 
maybe even more. The interstellar H2 and H-Lyman lines effectively remove a large fraction of the 
stellar flux shortward of 1020 A. Despite the plethora of interstellar lines in the FUSE spectrum of 
this star, we focus in this work on the FUV stellar wind features. Note that we did not attempt 
a detailed modeling of the FUSE spectrum. We did consider a model of the FUV spectrum of 
WR 46 (Crowther, private communication), but the direct comparison with the FUSE spectrum 
proved to be difficult given the wealth of interstellar absorption lines, so we did not try to fine-tune 
this model. It is however worth mentioning tha t the slope of the FUV continuum in the reddened 
model (E(B-V)=0.34, see lCrowther et al.lll995l ) was in rough agreement with the FUSE spectrum. 



WR 46 (WN3p) is the earliest WN subtype covered in the FUSE atlas of WR stars of 
Willis et al.l (|2004l ). The substantial interstellar contamination makes it very difficult to identify 
the expected S VI AA933.4, 944.5, C III A977, and N III A991 resonance lines. Iwillis et al.l J200J 
noted that the latter two are clearly not strong in emission or as P Cygni profiles, but that the S VI 
doublet seems to be present as a P Cygni profile. We can indeed identify the absorption troughs 
of this P Cygni profile around ~925 A and ~935 A. Note that this S VI resonance doublet is also 
observed as a P Cygni p r ofile i n the FUSE spectrum of the other WN3 star (HD 32109, WN3b) of 
the atlas of I Willis et aD (|200J). HD 32109 is in the LMC and suffers far less ISM contamination, 
which makes the identification of this spectral feature much easier. 



Willis et al.1 (|200J) also reported that the He II A1085, P V AA1118, 1128, and possibly Fe VI 



A1168 lines appear present as weak emissions or P Cygni profiles in WR 46. We also note the 
presence of a relatively wide absorption trough around 1060 A which might be the absorption 
component of a currently unidentified weak P Cygni profile. The FUSE spectrum of WR 46 
is dominated by a very strong saturated P Cygni profile in the O VI AA10 32, 1038 res o nance 
doublet. Given the high level of ionization potential of species in the wind, IWillis et al.1 ([2004) 
attributed the formation of O VI in this case to normal photoionization as opposed to Auger- 
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ionization formation in shocked wind gas. More generally, they considered the "super ions" of O VI 
and S VI as photoionized wind features for WN3-WN6 stars, while they are probably the result of 
Auger ionization in WN7-WN9 stars, and probably absent in WN10-WN11 stars. 

The presence of species of high ionization potential (138 eV for O VI) in the FUSE spectrum of 
WR 46 is not surprising given the very high effective temperatu re of this star. It is also in agreement 
with what is observed in its IUE spectrum \ Willis et al. 1986 ) which is dominated by N V AA1238, 
1242 and O V A1371 P Cygni profiles and He II A1640 emission, but in w hich no significan t N II I 
emission is detected. IUE observations of WR 46 were also presented by ICrowther et al.l (|1995l ) 
who noted that the He II A1640, N IV A1718, and C IV AA1548, 1550 lines are observed as pure 
emission features, with N IV and C IV extremely weak, while lines from higher excitation stages 
such as N V AA1238, 1242 and O V A1371 are observed as P Cygni profiles. The optical spectrum 
of WR 46 is also do minated by emission lines from the high-ionization species N V and O VI, and 
also from He II (e.g. lOliveira et al.l 12004 ). 



4.2. Time-series analysis of FUSE spectra 
To statistically determine which parts of the FUSE spectrum of WR 46 sh ow significant spectra l 



line variability, we performed a temporal variance spectrum analysis (TVS; iFullerton et al.l 11996) 
on the spectra from all the different observations summarized in Table [2j We did this separately 
for each FUSE channel, which allowed us, by comparing the results, to distinguish more easily 
real variable features from spurious variations. It was also easier to assess the quality of individual 
spectra by proceeding channel by channel. The co-addition of several FUSE spectral segments (with 
two to four segments overlapping depending on the wavelength region) introduces more variations 
of the continuum noise across the spectrum, which would complicate the determination of the 
weights in the calculation of the TVS (see below). Prior to computing the TVS of each channel, 
w e rebinned the spectra onto a common linear wavelength scale with a 0.13 A bin size. We refer 
to IFullerton et al.l (|1996l ) for the details of how to compute the TVS. We simply note here that the 
continuum noise was estimated by fitting a linear function to selected regions of the continuum for 
each channel, and that the errors on the flux of each pixel were taken from the flux error estimates 
output by CalFUSE, obviously modified according to the rebinning of the data. 

We present in Figure [2] (bottom panels) the TVS 1 / 2 spectrum for each channel. The quantity 
TVS 1 / 2 scales linearly with the size of the spectral flux deviations and is therefore more practical to 
use than the TVS itself because it gives a direct estimate of the amplitude of the spectral variations. 
Note that when plotting the TVS, we masked the wavelength intervals where peaks were caused 
by variability in airglow emission lines. We also masked spurious peaks caused by isolated outlying 
data points. These were common in the wings of narrow interstellar lines, since the TVS is very 
sensitive to small residual wavelength misalignment in the presence of steep spectral gradients. 

The dashed horizontal line plotted with the TVS 1 / 2 spectrum of each channel indicates the 
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99% statistical significance value for that channel. We can see from Figure[2]that there is significant 
variability in the blue edge of the O VI AA1032, 1038 P Cygni absorption components and also, to 
a lower level, in the S VI AA933, 944 P Cygni absorption components. These lines are in fact the 
two strongest P Cygni profiles present in the FUSE spectrum of WR 46. To show the evolution 
of these variable features with time for all FUSE observations of WR 46, we present in Figure [3] 
a gray-scale plot of the differences between each individual spectrum (with the spectral segments 
from all channels combined) and the reference time-averaged spectrum for the wavelength regions 
where significant variability was detected. In this figure, the bottom panels display each difference 
profile as an intensity profile. The top panel presents the time-averaged spectrum (thick line) with 
the spectra from all exposures overplotted (thin lines). For each line, we also overlay a velocity 
scale to show the velocity of the variable feature with respect to the center of the line. 

We can see that the variability in each component of the S VI doublet P Cygni absorption 
occurs mainly for velocities ranging roughly from -0.6 to a little faster than -v^, while the 
variability in the saturated O VI P Cygni absorption occurs in excess of -■u 0O , up to approximately 
-1.25 Voq. From the gray-scale plot of the consecutive observations of March 2006, we also see 
that the spectral features evolve on time scales of several hours, although the variability pattern 
is clearly not strictly periodic. It is indeed difficult to conclude anything about this variability 
pattern just by looking at the sequence of spectra. We do not clearly see displacement of spectral 
features in velocity space as a function of time, especially for the S VI doublet. The variations 
in each component of the S VI doublet absorption seem more or less correlated with each other. 
When the variable features are different in the two components, it is possibly due to the fact 
that the velocity intervals severely affected by interstellar absorption are also different for the two 
components. Finally, by comparing the evolution of the O VI and S VI profiles, we see no obvious 
correlation between the variability patterns in these features. 

Recall that for the above analysis, all the spectra were scaled to a common continuum level 
to focus on spectral-line variability. We can see from this analysis that the spectral lines do not 
seem to vary in strength over the years spanned by our FUSE observations. To identify possible 
long-term variability in the continuum, we looked at series of spectra that were not scaled to a 
common continuum level and found no evidence for significant long-term variations in the FUV 
continuum. 

4.3. FUSE light-curves and period search 

The FUSE light-curves extracted with the method described in <|3] are presented in Figure H] 
and Figure [5j Times are given in hours from the beginning of each observation (from the beginning 
of observation E1130103 for the March 2006 observations). As previously mentioned, these light- 
curves were extracted to provide a measure of the FUV continuum variability and of the integrated 
variability in selected spectral regions. Since the spectral variations generally appeared as broad 
features in the gray-scale plot of Figure [3], we believe that the light-curves for which the count rate 
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across the whole variable spectral region is integrated represent a good measure of the global line 
variability. 

As we can see from Figures H] and the amplitude of the FUV continuum light-curve is 
typically around ±5%. The amplitude of the O VI P Cygni absorption blue edge light-curve is 
typically around ±20%. The situation of the S VI P Cygni absorption light-curve is similar to that 
of the O VI light-curve, although the velocity intervals over which the changes occur are different. 
Interestingly, the amplitude of the variability in the P Cygni absorption components for both the 
O VI and S VI lines seems larger when the amplitude of the FUV continuum variability is also 
larger. We note some similarity in the shape of the three light-curves, although they do not appear 
strictly correlated. We will explore possible correlations between the different light-curves in £14.61 



To sear ch for periodicity i n the variability pattern of WR 46, we applied the fast Lomb-Scargle 



algorithm of lPress Sz Rvbickil (|1989l ) to the light-curves of the consecutive March 2006 observations. 
The resulting periodograms are displayed in Figure [6j We also plot on each periodogram a dashed 
line corresponding to the 99% significance level for variability based on white noise simulations. 
The FUV continuum periodogram shows a significant peak at a frequency of 3.15 cycles day^ 1 
(P=7.6±0.4 h), and another minor peak barely exceeding the 99% significance threshold at a 
frequency of 2.30 d _1 (P=10.4±1.4 h). The O VI P Cygni blue edge periodogram shows a significant 
peak at a frequency of 2.93 d^ 1 (P=8.2±0.5 h), and a second significant peak, just slightly weaker, 
at a frequency of 1.55 d _1 (P=15.5±2.5 h). This second peak might be linked to the first one since 
its frequency is almost 1/2 of the frequency of the main peak (the slight difference is perhaps just 
an effect of the limited time coverage of the data) . The situation is less clear for the S VI P Cygni 
absorption periodogram, partly due to the poor time coverage of the S VI light-curve for observation 
E1130105 (during which SiClB and SiC2A greatly suffered from pointing losses). There are two 
peaks around a frequency of 3 d _1 (8.00 h), but they just only slightly exceed the 99% significance 
level. 

The errors on the above periods were estimated as the HWHM of gaussian fits to the peaks 
in the Lomb-Scargle periodograms. Given the large widths of these peaks and the modest size of 
the data set from which we computed the periodograms, we will refer to the above measurements 
as dominant time scales and not as periods. From the results above, the variations in the FUV 
continuum and in the S VI and O VI P Cygni profiles are all consistent with a dominant time scale 
of ~8 h. 



4.4. XMM-Newton spectrum 

Even with the modest energy resolution of the EPIC instruments, it is clear from Figure [7] that 
the X-ray spectrum of WR 46 is dominated by emission lines. Given the weak N IV and strong N V 
lines in the optical and UV, the presence of N VI and N VII is not unexpected and it is clear that 
these ions account for most of the peak near 500 eV. N VI and N VII appear much stronger than 
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the neighboring O VII and O VIII lines near 600 eV, where the spectrum reaches a local minimum. 
The X-ray O/N ratio, althou gh uncertain becaus e of the poor resolution, looks roughly consistent 
with that seen in the optical. iGosset et al.l (|201ll ) modeled the X-ray spectrum with a combination 
of three optically-thin variable-abundance equilibrium plasmas of different temperatures. We took 
the alternative approach of building a model from individual ions from the K-shell emission of 
H-like and He-like ions and L-shell emission from Fe XVII- XVIII. As shown in Figure 6, this 
reproduces the observed spectrum well, giving the best-fit line intensities shown in Table El This 
table presents one entry per ion for the strongest line. There are many more lines in the model, 
but EPIC cannot resolve them. For Fe XVII and Fe XVIII, which both contribute many lines to 
WR 46, the relative intensities of all the lines were fixed to that of the strongest line, giving one free 
parameter. Similarly, EPIC cannot resolve the He-like triplets, so there is also one free parameter 
only each for N VI, O VII, Ne IX, Mg XI, Si XIII, S XV, Ar XVII. Because of the resolution, the line 
fluxes are only intended to be a rough guide to what future instruments might find. 



4.5. XMM-Newton light-curves and period search 



Figure [8] displays the simultaneous X-ray and UV light-curves obtained in the present work. 
We can see that the amplitude of the X-ray light-curve is around ±20%. The UV light-curve has 
quite an irregular shape and shows variations o n the order of ±5%. We again applied the fast 
Lomb-Scargle algorithm (jPress fe Rvbickil Il989l ) to search for periodicities in these light-curves. 
The resulting periodograms are displayed in Figure [91 with dashed lines still corresponding to the 
99% significance level for variability. The UV periodogram shows a significant peak at a frequency 
of 2.84 d _1 (P=8.5±1.8 h), while the X-ray periodogram shows a significant peak at a frequency 
of 2.66 d" 1 (P=9.0±2.5 h). 



From the XMM-Newton data, IGosset et al.l (|201ll ) extracted light-curves in various energy 
bands. They found a significant period of 7.9 h in the very soft X-rays (0.2-0.5 keV) and no 
significant variability at higher energies (above 0.5 keV). Our ti me scale for the X- ray variability 
in the total XMM-Newton band is consistent with that found bv lGosset et al.l (|201ll ) for their 0.2- 
0.5 keV light-curve, not surprisingly since the very soft photons contribute to a large proportion of 
the counts in our EPIC-pn light-curve. Since the X-ray spectrum is dominated by lines ( §4.4j) . the 
variability seen in the soft X-rays is also in lines, although it is not really possible to tell whether 
this is due to emission or absorption and which lines in particular are variable. 

Given the small number of cycles covered by the observations and the large widths of the 
peaks of the periodogram, we can only conclude that both the X-ray and UV light-curves vary on 
a dominant time scale close to 8-9 hours. These results are not incompatible with the variability 
time scale observed in the FUV continuum and P Cygni absorption components. We stress however 
that the two data sets are separated by more than four years. 
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4.6. Cross-correlation analysis 

Using the light-curves shown in the previous sections, we checked for possible correlations and 
delays between the FUV continuum and the P Cygni absorption component variations, and between 
the X-ray and UV variations. Such an analysis could eventually be useful to constrain a detailed 
model of the variability of WR 46. To test the overall similarity of the light-curves as a function of 
time lag, we shifted in time the simultaneous light-curves with respect to each other and calculated 



the correlation coefficient r (jPress et al.lll986l ) using the IDL routine C_CORRELATE.PRO. Before 
we did so, we rebinned the XMM-Newton X-ray and UV light-curves into common time bins. All 
the FUSE light-curves already had common time bins. 

In Figure [TU1 the cross-correlation functions are shown for the different light-curves compared. 
In the first case, we shifted the X-ray light-curve by At and calculated the correlation with the OM 
UV light-curve. In the other two cases, we shifted the O VI and the S VI light-curves by At and 
calculate the correlation with the FUV continuum light-curve. Note that this means that delays 
with respect to the UV or FUV variability correspond to negative values of At is this figure. 

We see from the plots of Figure [10] that consecutive positive correlation peaks are roughly 
spaced by the dominant variability time scale of 8 hours, and the same thing is true for consecutive 
negative correlation (anti-correlation) peaks. The cross-correlation analysis also tells us that the 
O VI, S VI and X-ray light-curves are more similar in shape to the FUV/UV light-curve when 
a shift is applied, suggesting that if they are related, there are possibly time delays between the 
variations in the continuum and the variations in the line- formation regions. 



5. Discussion 



5.1. FUV spectroscopic variability and CIRs 

The variability seen in excess of Vqo in P Cygni profiles is usually attributed to variable amounts 
of rarefied gas accelerated to velocities higher than . The shocks associated to this velocity excess 
can either originate from s mall-scale stochastic fluctuations coupled t o the radiative instability 



intrinsic to line-driven winds (Qwocki et al. 



1988 



Gavlev k, Owockilll995l ) , or fr om large-scale quasi- 



perio dic variability that can be induced by changes in the underlying star (e.g. ICranmer fc Qwocki 
199fih . 



The variability time scale of ~8 h that we found in the blue edge of the O VI doublet (and 
the S VI doublet to some extent) in WR 46 is close to the photometric and spectroscopic periods 
previously reported. It is also consistent with the variability time scale observed in the FUV 
continuum, suggesting that the FUV spectroscopic variability is not stochastic but caused by a 
large-scale perturbation controlled by an underlying quasi-periodic clock. 



CIRs are relatively large-scale spiral-shaped density and velocity perturbations that extend 
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from the stellar surface to possibly several tens of stellar radii into the stellar wind (|Mullanlll984l ). 
The Discrete Absorption Components (DACs) observed in the (unsaturated) absorption part of 
the UV P Cygni profiles in hot stars are commonly assumed to be a manifestation of CIRs. If the 
short-term variability of WR 46 is associated to CIRs, we might expect to see Discrete Absorption 
Components (DACs) in the unsaturated S VI P Cygni profile of WR 46. However, given the low 
signal-to-noise, sparse time coverage, and wealth of strong interstellar lines in this spectral region, 
it is not possible to identify in the S VI doublet the velocity displacement of an absorption feature 
which is typical of DACs (see Fig. [3]) . 

In addition to producing DACs, modeled CIRs can also modulate t he blue edge of the absor 
tion component of saturated P Cygni profiles on the same time scale ( Cranmer Owockilll99' 
This is in line with observations showing that th e steep blue edge of saturated UV lines varies 
in all cases when DACs are found in other lines (IKaper et al.l ll996). Based on the modulations 
observed in the blue edge of the O VI P Cygni absorption component of WR 46, we suggest that 
the short-term variability of the star might be related to CIRs, even though a typical DAC behavior 
could not be clearly identified in the unsaturated S VI line given our limited dataset. It is worth 
noting the case of the WR star EZ CMa (WR 6), for which 16 consecutive days of monitoring with 



IUE (>4 cycles) helped to link the observed UV P Cygni profile variability to CIRs (jSt-Louis et al 



19951 ). For this star, the variability is mainly found in the blue edge of the P Cygni absorption 
troughs. 



5.2. Variability originating from the photosphere 

The recent model atmosphere of WR 46 by Crowther (private communication) shows that the 
optical depth of the stellar continuum at 1000 A reaches a value of ~1 at the very base of the wind, 
and similarly for the 5000 A continuum. The FUV and optical continuum photometry, which both 
show variability at a ±5% level, are therefore dominated by the deep layers of the stellar wind, 
close to the hydr ostatic core . The a bsence of core eclipses in the FUV ( §4.3j) and optical continuum 



light-curves (e.g. I Veen et al.ll2002al ) suggests that the variability itself also originates from the very 



base of the wind and not from a binary companion. 

The similar periods or time scales (~8 h) found here or in previous studies in the optical, UV, 
FUV, and X-ray photometric and/or spectroscopic variability suggest that these are all related in 
some way. Even though they are not in phase (albeit correlated), the fact that the amplitude of 
variability in the S VI and O VI FUV lines is larger when the amplitude of the FUV continuum 
light-curve is also larger (see Fig. [3j 0] and [5]) also suggests a link between the continuum variations 
and line variability further out in the wind. The variability of WR 46 could originate from the 
photosphere of the star and propagate out, modulating the absorption of X-rays or the structure 
of wind shocks on the same time scale as the ~8-h period. 



If this is correct, we might observe delays between the continuum variations and the other 
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manifestations of the variability. According to the CMFGEN model of WR 46, the peak of the line 
formation region of the O VI doubl et is around 2 .5 while the formation of the S VI doublet 
is much closer to the photosphere. iGosset et al.l (|201ll ) also estimated that the N VI and N VII 
X-ray lines, which we argued dominate the X-ray spectrum and the X-ray variability (see §4.4[) . 
form above a few stellar radii. Given the parameters adopted in Table [TJ we should expect the 
line variability (in particular the O VI and X-ray variability) to lag the continuum variability by a 
few hours if perturbations are advected by the wind. The fact that the cross-correlation functions 
for the O VI and X-ray light-curves (Figure I10p do not show an extremum when the time lag is 
zero suggests that these are not in phase with the continuum variability. The presence of extrema 
corresponding to delays of a few hours are also not inconsistent with the expected time lags. 

It is worth noting that iBerghoefer et al. ( 1996 ) found correlated variability in the X-ray and 
Ha emission from the 04If supergiant £ Puppis. They suggested this was evidence for periodic 
density variations propagating from the base of the wind to further out where X-rays are produced, 
an interpretation that is very similar to what we propose for WR 46. 



5.3. Lack of stability of the period 



We did not find any significant long-term changes in the FUSE spectra between 2000 and 
2006 (in continuum and line strength). As a consequence, we found no evidence for a change of 
the mass-loss rate , unlike the change in the global mass-loss rate in the early 90's reported by 
Veen et al.l (|2002bl ). The variability time scales of our 2 006 FUV data are c onsistent with the 
dominant photometric and spectroscopic periods found by lQliveira et al.l (120041 ) in their data from 
1998-1999. It is therefore possible that while the global mass-loss rate did not change, the dominant 
period of the system also stayed the same between 1998 and 2006. This is not incompatible with 
the suggestion t hat the long-term changes (luminosity, mass-loss rate) in WR 46 are linked to the 



period changes (IVeen et al 



2002 



3D). 



Inspection of Figure Q] reveals that there is no clear trend (ei ther increas i ng or d ecreasing) in 
the evolution of the period of WR 46. From two intense data sets, I Veen et al.1 (|2002al ) showed that 
there is a decrease of the photometric period from 6.78 h to 6.54 h between 1989 and 1991, and 
that this period change is significant. Then, most of the periods reported for later years are longer, 
although it is hard to assess the reliability or precision of some of these periods, especially when no 
uncertainties were quoted. Interestingly, the time scales that we measured from the FUV continuum 
and O VI blue e dge light-curves (see ^4.3f) are significantly longer than the dominant photometric 
periods found by I Veen et al.l (|2002al ). Even when considering the rather large uncertainties on our 
measured time scales, our values are not consistent with the periods from 1989 and 1991. 

This apparent lack of stability of the period of the system is a strong argument against rota- 
tional modulation as the single origin of the short-term variability of WR 46. A well-known example 
where rotational modulation has been linked to the short-term variability is the case of EZ CMa 
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(WR 6) OSt-Louis et al.lll995l : Morel et al.lll998l : iFlores et all 120071 ). This star always displays the 
same period of 3.766 days, even though its variability is epoch-dependent. If rotational modulation 
is responsible for the variability, the period should indeed always be the same or should change very 
little and gradually over time. We can also argue that rotation of a single WR star is not a likely 
cause of short-term variability with a period of a few hours because it would generally require a 
rotational velocity that is high relative to the break-up velocity of the star. For WR 46, assuming 
a rotation period of 8 h and the parameters of Tabled! we find that the star would rotate at about 
40% of its critical rotational velocity at the equator, with a velocity of ~ 440 km s _1 . This is not 
im possible per se but still rather h igh compared to the first estimates of rotation rates of WR stars 
by 



Chene fc St-Louid (|2008l . |2010| ) who found rotational velocities typically of a few tens of km s 



The lack of stability of the period is also difficult to explain in terms of binarity alone. After 
considering c onservation l aws (s pin and orbital angular momentum, energy) and various exchange 
mechanisms, IVeen et al.l (|2002a ) did not exclude that the observed period change from 1989 to 
1991 represents an orbital change due to the spiral-in of a companion, but they considered such 
a scenari o highly unlikely. Between 1989 and 1991, the decrease in the photo metric period is 1 4 
minutes (jVeen et al.ll2002al ). which is extremely large for any binary system (see lVeen et al.ll2002d ). 
This spiral-in scenario seems even more unlikely when we consider the subsequent measurements 
of longer (and possibly multiple) periods. Again, the time scales that we measured from our 2006 
FUSE obser vations are s ignifica ntly longer than the dominant photometric periods found in 1989 
and 1991 bv lVeen et aD (j2002al ). 



5.4. X-rays and the Hatchett-McCray effect 



If we compare the X-ray properties of WR 46 with those of the three WR+compact sys- 
tems known, it seems very unlikely that WR 46 belongs to the same class. The most famous 
WR+compact sys tem i s probably the strong galactic X-ray source Cyg-X3 (L v ~ 10 38 erg s~ 



Sch mutz et al.l e.g. ll996l . , and references therei n) . It has an orbital period of 4.8 h (jParsignault et al 



19721 ) and is t hought to consist of a WR star (Ivan Kerkwiik et al 



19921 ) possibly accompanied by 



a black hol e (jSchmutz et al.1 119961 ; Hanson et al.l l2000l ). although this has been contested (e.g. 
Mitral Il998l ). A firs t extragalactic WR+compact candidate, IC 10 X-l, was detected in the star- 
burst galaxy IC 10 (jBauer Brandtll2004l : IWang et al.ll2005l ). It also has a high X-ray luminosit y 



(Ly(0.2 - 10 keV) ~ 1.2 x 10 38 erg s _1 ), and its period is close to 35 h (jPrestwich et al J 120071 ). 



Silverman Sz Filippenkd (120081 ) confirmed that this syste m is a WR + bl a ck hole binary. Another 
extragalactic candidate, NGC 300 X-l, w as identified by Carpano et al. ( 2007al jbl) and then 



_ con- 
firmed as a WR + black hole binary by iCrowther et al. (2010). Its observed luminosity in the 
0.2-10 keV band is ~ 2 x 10 38 erg s~ 1 , and it was found to have a period close to 32 h. As well as 
the orders of magnitude higher luminosities and much harder spectra, those accreting X-ray sources 
show a much smoother continuum spectrum compared to the line-dominated X-ray spectrum of 
WR 46 (Figured). Again, we stress that the X-ray spectrum and luminosity of WR 46 are typical 
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of single WN stars. 



Estimating the total X-r ay luminosity produced by Bondi-Hoyle accretion of a stellar wind 



onto a degenerate object (e.g. IStevens Willis! Il988l ) for the parameters of WR 46 (assuming a 



circular orbit with P=8 h), we obtain Lx ~ 10 ergs s . This expected X-ray luminosity is several 
orders of magnitude higher than what is observed for WR 46, and this is very unlikely caused by 
wind attenuation of the X-rays. For example, the predicted X-ray accretion luminosity of EZ CMa 
is ~ 10 36 erg s _1 , about 3 orders of magnitude higher than what i s observed, an argu ment that 



has been used against an acc reting neutron star companion scenario ([Skinner et al.lll998l ). For this 



star, IStevens &: Willis! (|1988l ) performed numerical calculations accounting for wind attenuation for 



a range of system parameters and they predicted values of Lx ~ 10 — 10 erg s , still much 
larger than the observed value. WR 46 is a weak-lined WN star and its wind is less dense than 
that of EZ CMa, so the effect of wind attenuation should not be stronger than in EZ CMa. Thus, 
it is highly unlikely that the wind of WR 46 hides a very luminous X-ray companion, especially one 
that would imply a >4 orders of magnitude difference between the intrinsic and absorbed X-ray 
luminosity. 

V arious suggestions like centrifugal inhibition of accretion (llllarionov Sunvaevlll97.4 l Stella et al 



19861 ) have been discussed as pote ntial mechanisms t hat could explain the lack of X-rays from a 
compact companion. For WR 46, 1 Veen et al.l (12002a ) also discussed other possible types of com- 
panion like a low-mass main seq uence star, which would however imply an unexpected initial mass 
ratio (e.g. iGarmany et al.lll980l ). Some of these more unusual types of companions might not be 
very luminous in X-rays. However, in such cases, even though the observed level of X-rays could 
possibly be reconciled with the presence of a companion, its impact on the stellar wind would 
remain minimal as we illustrate below. 

We have already reviewed many problems associated with the interpretation of the variability of 
WR 46 in terms of binarity alone. Nevertheless, assuming that more than one cyclical mechanism 
could be at play in this star, it is worth checking if a companion could be responsible for the 
variability observed in the P Cygni profiles. An X-ray source immersed in the stellar wind of its 
companion star can ionize a surrounding region and lead to observable variations. In particular, such 
a system could show variability with orbital phase in the P Cygni profile of ultr aviolet resonance 
lines, an effect known as the Hatchett-McCray effect (jHatchett Sz McCravl 119771). w hich has been 



observed in several high-mass X-ray binaries (HMXRBs) (e.g. 



van Loon et al 



2001 



If an X-ray source was present in the wind of WR 46, then there would be a region around 
this source in which atoms are ionized to higher stages. In this region, the fraction of O VI or 
S VI with respect to other ionization stages would change. For a given transition, an increased 
number of scatterers in the ionization region would lead to an enhanced P Cygni absorption and/or 
emission, depending on where the ionization region is with respect to the WR star and the observer. 
Similarly, if the number of scatterers decreases in the region surrounding the X-ray source, the 
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P Cygni absorption and/or emission weakens. Note that it is possible for the variations in the 
P Cygni emission component to remain unnoticed while variations are observed in the absorption 
component, because the region ionized by the X-ray source represents a smaller fraction of the 
volume contributing to the emission component compared to the fraction of the volume contributing 
to the absorption component. 

Several factors however appear to go against an interpretation of the variability of the P Cygni 
profiles of WR 46 in terms of the Hatchett-McCray effect. First of all, the X-ray luminosity of 
WR 46 (~10 32 erg s _1 ) is very low compared to the typical X- ray luminosity of kno wn HMXRBs 
showing the Hatchett-McCray effect (~10 35 -10 38 erg s" 1 ) (e.g. Ivan Loon et alJboOll ). Also, since 



the size of the ionized region for a given species and ionization stage sc ales linearly with the X-ra y 



luminosity of the companion and inversely with the density of the wind (jHatchett &: McCraylll977l ). 
we expect this region to be very small for WR 46 which has a relatively low X-ray luminosity and a 
denser wind t han typical O stars. We ca n estimate the size of such an ionized region by considering 



the quantity (jHatchett fc McCrayHl977l ) 



/ \ 2 

Lx 4-7T Lx m I r 



€(r,rx) = = -r- v(r) — , (5) 

n(r) r x M \ r x J 

where n(r) is the local number density of the gas, rx is the distance from the X-ray source, m is the 
average mass per ion (we assume here a pure helium atmosphere), Lx is the intrinsic (unabsorbed) 
X-ray luminosity of the companion, v(r) is the velocity of the WR wind, and M is the WR mass-loss 
rate. In the limit £ — > 0, the ionization balance of the material is unaffected by the presence of 
the X-ray companion and t hus mainly governed by the radiation field of the WR star. According 



to the ionization models of iKallman fc McCravl (|1982| ). significant changes in ionization fraction 



generally only occur for values of log£ > 1.6. The variations of log £(r, rx) in the orbital plane for 
the case of WR 46 are shown in Figure [TT] for an assumed intrinsic X-ray luminosity of 10 33 , 10 34 , 
or 10 35 erg s . For this calculation, we assumed an orbital period of 8 h and a circular orbit, a 
/3 = 1 velocity law, and we used the values of Table [TJ The corresponding distance of a compact 
companion is 5.3 Rq, or 1.83 R*. 

The Hatchett-McCray effect is obviously more difficult to observe if the ionized region is 
small. If it is small enough, variability should be limited to the absorption components of P Cygni 
profiles and to the fraction of the orbit during which the companion is between the star and the 
observer. We can see from Figure [TT] that even in the rather unlikely case where the unabsorbed 
X-ray luminosity of the companion is 10 35 erg s _1 (see discussion above), the region within which 
log^ > 1.6 is relatively small, and it is unlikely to perturb the column in front of the star for most 
of the orbit. Note also that equation [5] is valid for an optically thin gas, and that in reality the mean 
path length of the X-ray photons should be considerably shorter. The X-ray ionized zone is thus 
likely even less extended than what is shown in Figure [TT] The fact that we see spectral variations 
during a major fraction of the ~8-hour cycle appears in contradiction with the Hatchett-McCray 
effect of a small ionized region. The ionized region would also have to be much larger than predicted 
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here in order to encompass the region of the wind between -0.6 Voo and -i>oo> where variability is 
observed in the S VI P Cygni absorption. 



5.5. Nonradial pulsations 



The variability observed in this work in the FUV continuum and in the UV, as well as the 
optical continuum variability seen by other authors, would find a natural explanation in the context 
of pulsations. It could either be caused directly by brightness variations at the stellar surface, or 
by velocity /density perturbations at the stellar surface deforming the continuum forming layer. As 
discussed in §5.21 this could then induce variability further out in the wind. 

The scenario of non-radial p ulsations is also appealing to explain the multiple photometric 
periods present at a given epoch (|Veen et alj|2002al : lOliveira et al.ll2004l ). Moreover, the different 
periods measured and the fact that the variations are not strictly periodic might easily be explained 
by the presence of multiple pulsation periods. The interference patterns between different modes 
could result in a complex light-curve, varying both in shape and amplitude. It is also conceivable 
that a time-limited data set would allow to identify only one of the periods, and that another 
time-limited data set would show another of the periods. 

Many studies have explored and discussed the possibility of pulsations in Wolf-Rayet stars 



Ce.g.lNoels Scuf lairc 



Dorfi et al 



1986; 



2006; iGlatzel 



Glatzel et al. 



200. 



1993 



For WR 46, Veen et al 



Lefevre et al 



2005 



Townsend fc MacDonald 



2006 



(j2002d ) proposed a simple geometrical 



model in which non-radial pulsations would be able to cause radial-velocity shifts of the wind 
emission lines. They suggested that the lowest order sectoral mode 1=1 and |m|=l, which results 
in a single bright hemisphere travelling along the equator, could produce a distorsion of the wind 
of WR 46 (see their Fig. 5) and lead to a double-wave photometric period and a single-wave 
radial-velocity period, as observed. The bright hemisphere could distort the wind by causing an 
asymmetric density flow, with the "one-armed" wind appearing to rotate as the pulsational pattern 
travels over the stellar surface. In such a scenario, non-radial pulsations could generate CIRs in 
the wind. 



6. Conclusion 

From our time-resolved FUSE observations of WR 46, we detected "blue edge" variability 
in the O VI doublet P Cygni absorption and line-profile variations in the S VI doublet P Cygni 
absorption. Both variations occur on a time scale of about 8 h, consistent with the time scale 
found in the FUV continuum light-curve. A similar time scale was also recovered in UV and X-ray 
light-curves from archival XMM-Newton observations. 

We noted that the X-ray spectrum and luminosity of WR 46 are typical of a single WN star. 
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If a companion with a lower X-ray luminosity like a rapidly rotating neutron star was present, its 
effect on the stellar wind would be minimal, and its limited ionization effect (the Hatchett-McCray 
effect) could not explain the variability observed in the FUV P Cygni profiles. 

We also pointed out that the fact that the observed period is not stable (changing with time 
with no particular trend) is a strong argument against rotational modulation or binarity alone as 
the origin of the short-term variability. 

Circumstantial evidence suggests that the variability of WR 46 is rooted in photospheric per- 
turbations, and that CIR-like structures might induce the observed P Cygni line profile and X-ray 
line variability. We proposed that non-radial pulsations might be the source of these perturbations, 
although an improved theoretical understanding of non-radial pulsations in WR stars and how they 
propagate in the stellar wind is required to really investigate this possibility further. 

Meanwhile, we suggest that additional intensive and long-term monitoring campaigns of WR 46 
would certainly prove useful, for example to follow the evolution of the period or the recurrence of 
multiple periods. If possible, continuous photometry for several days or weeks would be a perfect 
way to test the scenario of non-radial pulsations with multiple periods. 

We are thankful to Paul Crowther for providing results from his recent wind model of WR 46. 
We also thank the anonymous referee for comments which helped to improve the paper. V. H.- 
B. acknowledges support from NSERC (Canada) for a postgraduate scholarship. N. S.-L. is also 
grateful to NSERC for financial assistance. 
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Fig. 1. — Periods reported for WR 46 plotted versus year of observation. The photometric periods 
are shown with stars and the spectroscopic periods with circles. The upper panel shows the domi- 
nant periods for a given data set and the bottom panel shows periods with weaker amplitudes. The 
references below are identifi ed on each measure ment by their numb e r. Er ror b ars are shown whe n 
ava ilable. References - (1) IVeen et al.1 d2002ah . f2\ loiiveira et al.1 d2004h. (3 ) IVeen et al.1 \ 19951 ). 



(4) iNiemela et al.l (|1995T ). (5) Marchenko et al.l (|2000|), (6) iGosset et al.1 (|201lh . (7) This work. 



-26- 



7.0 
6.0 
5.0 
4.0 
3.0 
2.0 

1.0 
0.0 

1.5 
1.0 
0.5 
0.0 



£ 

u 
a; 



a; 



SiClB 

: A 


average spectrum 


Tvg l/! 






"l 1 1 1 





920 



940 



960 



960 



7.0 
6.0 
5.0 
4.0 
3.0 
2.0 
1.0 
0.0 
1.5 

1.0 
0.5 

o.o t 

980 
7.0 
6.0 
5.0 
4.0 
3.0 
2.0 

1.0 
0.0 
1.5 

1.0 
0.5 
0.0 

7.0 
6.0 
5.0 
4.0 
3.0 
2.0 

1.0 
0.0 

1.5 

1.0 
0.5 
0.0 



LiF2B average spectrum i 
VI 



SiCIA average spectrum 
VI 




I | , , , 1 , , 1 1 , 1 1 1 r 

1020 1040 1060 1080 



L 


FIB aver 

ft 


age sped 
1 


rum 

A 




TVS 1/2 





i — 1 — 1 — 1 — r 

1100 1120 



-| — i — i — i — | — i — i — i — r 
1140 1160 1180 



1 

1 

1.0 - 



SiC2A average spectrum 

S VI 
I 1 



I 



TVS ' 



920 



1 

940 



1 

960 



9B0 



I 

1000 




LiFIA average spectijum 
VI 





SiC2B average spectrum 




" | TVS 1/2 




Wavelength (A) 



Fig. 2. — Top panels: Time-average spectrum for each FUSE channel. Bottom panels: Square- 
root of the temporal variance spectrum (TVS 1 / 2 ) of our time series of spectra. The dashed line 
indicates the 99% confidence level for variability. The wavelength intervals with significant spectral 
variability and for which light-curves were computed are highlighted in gray. 
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Fig. 3. — Gray-scale plot of the differences between individual spectra and the time-averaged 
spectrum (top, thick line) for the S VI AA933, 944 P Cygni absorption trough (left) and the O VI 
AA1032, 1038 P Cygni absorption blue edge (right). Time is given in hours since the beginning of 
each observation. The velocity scales above the top panel are given with respect to the center of 
the S VI A933, S VI A944, and O VI A1032 lines, using the adopted = 2775 km s _1 . The narrow 
absorption features are of interstellar origin. Smoothing factors of 15 and 5 were applied to the 
S VI and O VI difference spectra respectively. The real signal-to-noise of the difference spectra is 
therefore lower than seen here. Fluxes are in units of 10~ 12 ergs s _1 cm -2 A -1 . 
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Fig. 4. — Light curves for the consecutive March 2006 observations in spectral windows correspond- 
ing to the FUV continuum (top panels) and the variable regions of the O VI and S VI P Cygni 
profiles (bottom two panels). 
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Fig. 5. — Same as Figure0]but for the non-consecutive observations P1170601, Z9900101, E1130101, 
and El 130102 respectively from left to right. 
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Fig. 6. — Lomb-Scargle periodograms for the combined light-curves of March 2006 FUSE obser- 
vations, in spectral windows corresponding to the FUV continuum and the variable absorption 
component of the O VI and S VI P Cygni profiles. The dashed lines indicate the 99% significance 
level. 
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Fig. 7.— The XMM-Newton X-ray spectrum of WR46 with EPIC-pn in black and EPIC-MOS1 
and EPIC-MOS2 added together in gray. Data are shown with error bars while the continuous 
lines show a pure emission-line model from H-like and He-like ions of C, N, O, Ne, Mg, Si, S and 
Ar along with L-shell emission from FeXVII and FeXVIII. 
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Fig. 8.— Top panel: Broadband (0.3-10 keV) XMM-EPIC X-ray light-curve of WR 46. Bottom 
panel: Simultaneous XMM-OM UV light-curve using the UVM2 filter for the first half of the 
observation and the UVW2 for the second half. Note that before normalizing the UV light-curve, 
the measured rate in the UVM2 filter was converted to its equivalent in the UVW2 filter using the 
method described in <33.2l 
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Fig. 9. — Lomb-Scargle periodograms for the XMM X-ray and UV light-curves of WR 46. The 
dashed lines indicate the 99% significance level. 
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Fig. 10. — Linear correlation coefficient r as a function of the time translation applied to the light- 
curves. The top panel shows the correlation of the shifted XMM X-ray light-curve with the XMM 
UV light-curve. The middle and bottom panels show respectively the correlation of the shifted 
O VI and S VI FUSE light-curves with the FUSE FUV continuum light-curve. 
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Fig. 11. — Contours of constant ionization parameter log£ (in erg cm s _1 ) for an ionizing com- 
panion located at (0, 5.3) and having an intrinsic (unabsorbed) X-ray luminosity of 10 33 , 10 34 , or 
10 35 erg s _1 (from left to right). The Wolf-Rayet star is placed at the origin of the plots. The 
region within which log£ > 1.6 is indicated in light gray. 



-35- 



Table 1. Adopted stellar and wind parameters of WR 46 

Parameter Value 

R* 2.9 Rq 

T, 90,210 K 

M 18 M© 

L 0.5xl0 6 L Q 

Voa 2775 km s _1 

M 4xl0~ 6 M Q yr" 1 

d 4.0 kpc 

L x (0. 2-10.0 keV) 7.7 x 10 32 erg s _1 



Table 2. Log of FUSE observations 



Data set* Start time (UT) End time (UT) Aperture MODE Number of exposures t eX p(s) 



PI 170601 2000-04-08 03:36:26 2000-04-08 07:25:06 LWRS TTAG 4 4891 

Z9900101 2002-06-28 14:31:54 2002-06-28 15:56:39 LWRS TTAG 2 3946 

E1130101 2005-07-17 04:09:40 2005-07-17 10:26:59 LWRS TTAG 9 10210 

E1130102 2005-07-20 03:46:53 2005-07-20 08:12:58 LWRS TTAG 6 6199 

E1130103 2006-03-09 14:28:48 2006-03-09 20:23:59 LWRS TTAG 4 9955 

E1130104 2006-03-10 02:46:20 2006-03-10 19:38:25 LWRS TTAG 11 29682 

E1130105 2006-03-11 03:44:02 2006-03-11 17:06:18 LWRS TTAG 9 28213 



"The data are archived in the Multimission Archive of the Space Telescope Science Institute (MAST). 
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Table 3. List of X-ray lines and measured fluxes 



Ion Wavelength (A) Flux (cm 2 s 1 ) 



C VI 


33.7342 


1.8 


± 


1.0 


X 


10" 


-4 


N VI 


28.7871 


1.8 


± 


0.5 


X 


10" 


-4 


N VII 


24.7793 


7.2 


± 


0.4 


X 


10" 


-5 


O VII 


21.6020 


2.4 


± 


0.4 


X 


10" 


-5 


O VIII 


18.9671 


5.6 


± 


1.0 


X 


10" 


-6 


Fe XVII 


15.0150 


6.9 


± 


1.0 


X 


10" 


-6 


Fe XVIII 


14.2080 


2.4 


± 


0.6 


X 


10" 


-6 


Ne IX 


13.4471 


1.2 


± 


0.1 


X 


10" 


-5 


NeX 


12.1321 


2.7 


± 


0.4 


X 


10" 


-6 


Mg XI 


9.1688 


3.3 


± 


0.1 


X 


10" 


-6 


Mg XII 


8.4192 


1.2 


± 


0.1 


X 


10" 


-6 


Si XIII 


6.6480 


2.4 


± 


0.1 


X 


10" 


-6 


Si XIV 


6.1804 


9.6 


± 


2.1 


X 


10" 


-7 


sxv 


5.0387 


3.4 


± 


1.2 


X 


10" 


-7 


SXVI 


4.7274 


9.3 


± 


5.8 


X 


10" 


-7 


Ar XVII 


3.9488 


2.9 


± 


1.3 


X 


10" 


-7 


Ar XVIII 


3.7311 


3.6 


± 


1.0 


X 


10" 


-7 



